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Abstract. We discuss new high-field magnetization data recently obtained by Tsirlin et al. 
for layered vanadium phosphates in the framework of the square-lattice model. Our predictions 
for the saturation fields compare exceptionally well to the experimental findings, and the strong 
bending of the curves below saturation agrees very well with the experimental field dependence. 
Furthermore we discuss the remarkably good agreement of the frustrated Heisenberg model on 
the square lattice in spite of the fact that the compounds described with this model actually 
have a lower crystallographic symmetry. We present results from our calculations on the 
thermodynamics of the model on the orthorhombic (i.e., rectangular) lattice, in particular the 
temperature dependence of the magnetic susceptibility. This analysis also sheds light on the 
discussion of magnetic frustration and anisotropy of a class of iron pnictide parent compounds, 
where several alternative suggestions for the magnetic exchange models were proposed. 



1. Introduction 

The frustrated 5 = 1/2 Heisenberg model on the square lattice, the J1-J2 model, appears to 
describe well the thermodynamic and magnetic properties of two classes of vanadium compounds 



of type Li2VOX04 {X = Si, Ge) [T] and AA'VO(P04)2 {A, A' = Pb, Zn, Sr, Ba) [21 [3]. They 



consist of V-oxide pyramid layers containing V^"^ ions with S = 1/2. Prom the analysis of the 
temperature dependence of the heat capacity and the magnetic susceptibility in zero (or small) 
fields, the frustration ratio J2/J1 can be ontained However, an ambiguity remains with 
respect to the relative sign of the two exchange constants, which can be resolved by analyzing 



the behavior of these materials in finite fields [5] . The average exchange constants Jc = y + Jf 
of these materials are low enough such that their saturation fields are experimentally accessible. 
In this article we discuss the high- field magnetization of the J1-J2 model and compare our 
findings to recent measurements [6]. 

Although we describe the physics of the J1-J2 compounds using a square-lattice model, their 
true crystal structure corresponds to a two-dimensional lattice with lower symmetry. Therefore, 
we introduce an additional spatial anisotropy in the ab plane assuming orthorhombic symmetry, 

1. e., a rectangular lattice. 

2. Model Hamiltonian 

The effective Hamiltonian on the rectangular lattice has the form 
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Figure 1. Left: Comparison of saturation fields from ED and spin-wave theory with = ^, and 
from high-field experiments [6j. The compounds are (1) PbZnVO(P04)2, (2) Nai,5VOP04Fo.5, 
(3) Pb2VO(P04)2, (4) SrZnVO(P04)2, and (5) BaCdVO(P04)2- The agreement with the 
saturation fields for the columnar phase gives a direct proof that all compounds have CAF 
order. Right: Magnetization for BaCdVO(P04)2 Solid line denotes experimental data and 
filled symbols the data from T = Lanczos calculations for different cluster sizes. 



where Jia and Jib denote the nearest-neighbor exchange along the a and b directions, and J2 
labels the diagonal next-nearest neighbor exchange. A more convenient parametrization is 

Jia = \/2 Jc cos (j) cos 0, J2 = Jc sin 4>, 

Jib = V2Jc cos (/) sin 9, Jc = \j\ {Jia + Jfb) + Jl (2) 

introducing a frustration angle (p, an anisotropy parameter 9, and an overall energy scale Jc. 
For 9 = vr/4, the above Hamiltonian reduces to the square-lattice case. 

The results presented in the following paragraphs are obtained by applying linear spin- 
wave theory (LSW) and exact diagonalization (ED) for finite clusters at finite temperatures 
(FTLM) [ilSlE]. 

3. Saturation field 

From LSW the saturation field for the square lattice is given by the instability of the fully 
polarized state against a single-magnon excitation. This result is in exact agreement with ED 
for antiferromagnetic J2. However, exact diagonalization for our finite clusters reveals that for 
ferromagnetic Ji, a AS" = 2 instability determines the saturation field |5]. 

The left-hand side of Fig. [T] shows a comparison of the saturation fields determined by LSW 
and ED for the columnar (CAF) and Neel (NAF) antiferromagnetic phases with the experimental 
values determined from high- field measurements [6]- The predicted theoretical values are based 
on fits of our FTLM data and of a high-temperature series expansion [21 [3l El [8] to the 
temperature dependences of the low-field susceptibilities. The experiments agree surprisingly 
well with the predicted CAF values, demonstrating that all compounds order in a columnar 
magnetic structure at low temperatures. 

On the right-hand side of Fig. [H the field dependence of the magnetization for 
BaCdVO(P04)2 is displayed, together with zero-temperature data from our Lanczos calculations 
for different cluster sizes using a Bonner-Fisher construction [51 [H] • Given the small size of the 
clusters involved, the agreement is well, apart from low fields, where finite-size effects are most 
prominent. 




Figure 2. FTLM results for the magnetic susceptibility xi^)- The left-hand figures show 
the position, the right-hand ones the temperature where xC^) reaches its maximum. Top: 
Dependency on the frustration angle (j). The vertical lines distinguish the different classical 
phases, FM, NAF and CAF. Bottom: Dependency on the anisotropy parameter 9 at fixed 
frustration angle (p/ir = 0.625. The vertical line shows the isotropic case, 9 = 7r/4. 



4. Extension to the orthorhombic (rectangular) case 

Up to here, we discussed the layered vanadium phosphates in the context of the square- 
lattice Heisenberg model. However, their crystal structure corresponds to a lattice with lower 
symmetry [10] . We therefore investigate the impact of an additional anisotropy of the nearest- 
neighbor interactions in the ab plane, characterized by the angle 9 defined in Eq. [2l 

In Fig. [21 FTLM results for a tile of size 20 are shown. The position (left) and the value 
(right) of the broad maximum of the magnetic susceptibility xi^) are plotted as a function of 
the frustration {(f>) and the anisotropy (9) angles. The top curves correspond to the isotropic 
case (Jia = Jib, 9 = j), while the bottom ones show the effect of the anisotropy parameter 9 
for constant (p/ir = 0.625 (CAF regime). The change both in value and temperature is small. 
Therefore introducing an anisotropy within the columnar phase has comparatively little effect 
on the temperature dependence of xi^)- For further clarification, a contour plot of the ground 
state energy (right) of the Hamiltonian, Eq. [T] and the temperature of the maximum of x(^) 
(left) as a function of (j) and 9 are shown in Fig. [3l The model has four classical phases, one 
FM, one Neel AF, and two columnar AF phases along the crystallographic a and b directions. 
Inside the AF regions, the parameter dependence is weak. This explains the validity of the 
square-lattice Heisenberg model in describing the experimental results on the thermodynamics 
of the compounds, which have lower than tetragonal symmetry. 

It has been proposed that magnetic frustration is a key feature of the magnetic properties 
of ferropnictides. However, several different exchange models are discussed |1 H \12 \ [TH | I14j. in 
particular models with a spatial anisotropy as described here. As an example. Table [T] shows 
the experimental and theoretical values for moments and exchange constants in some iron 
pnictide parent compounds, which are all in the CAF regime. We conclude that the anisotropy 
of Jia,b stabilizes the CAF phase. In particular the values with Jih ~ and Jia/2J2 — 1 
(0 ~ 0, i;/) ~ 0.15-7r) correspond to the stable CAF region. It is obvious from Fig. [3] (left) that 




Figure 3. Left: Contour plot of ground state energy as function of anisotropy (9) and frustration 
{(p). Open circle designates usual NAF {6 = ^,0 = 0). Values in table 1 are represented by 
black (exp.) and white (theory) symbols. The white lines show the boundaries between the four 
classical phases, CAFa^^, NAF and FM. Right: Contour plot for the position of the maximum 
susceptibility. 
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0.19 
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0.22 
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0.11 


-0.04 


BaFe2As2 
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0.28 
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35 


39.1 


0.35 


0.25 


BaFe2As2 


m 


0.54 


36 


-7 


18 


31.6 


0.19 


-0.06 


CaFe2As2 


m 


0.75 


27.4 


-2.1 


14.5 


24.3 


0.20 


-0.02 


BaFe2As2 


m 


0.84 


36.1 


-2.6 


12.0 


38.0 


0.10 


-0.02 


SrFesAss 


m 


0.84 


35.3 


2.2 


13.4 


28.4 


0.16 


0.02 



Table 1. Fe pnictide moment /i = 2SfiB and exchange interactions (in meV) from experiment 
(top) and theory (bottom). 

these values are quite distant to the strongly frustrated point {6 = 0.257r, (j) = O.lSvr) where the 
CAFfl b and NAF phases meet. The moment reduction by quantum fluctuations for the former 
values is comparable to that of the simple unfrustrated NAF (open circle). 
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